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The quantum dynamics of a chaotic billiard with moving 
boundary is considered in this work. We found a shape pa- 
rameter Hamiltonian expansion which enables us to obtain 
the spectrum of the deformed billiard for deformations so 
large as the characteristic wave length. Then, for a specified 
time dependent shape variation, the quantum dynamics of a 
particle inside the billiard is integrated directly. In particular, 
the dispersion of the energy is studied in the Bunimovich sta- 
dium billiard with oscillating boundary. The results showed 
that the distribution of energy spreads diffusively for the first 
oscillations of the boundary ((A 2 E) — 2 D t). We studied 
the diffusion contant D as a function of the boundary veloc- 
ity and found differences with theoretical predictions based 
on random matrix theory. By extracting highly phase space 
localized structures from the spectrum, previous differences 
were reduced significantly. This fact provides the first numer- 
ical evidence of the influence of phase space localization on 
the quantum diffusion of a chaotic system. 

PACS numbers: 05.45. +b, 03.65.Sq, 03.20.+i 



I. INTRODUCTION 

The quantum dynamics of a classically chaotic system 
is in the present a topic of very active interest. Specially, 
with the great development in artificially fabricated small 
devices as quantum dots Q, where quantum manifesta- 
tion of classical chaos || plays an important role. 

Systems governed by a parameter dependent Hamilto- 
nian H(£) are excellent models to study quantum man- 
ifestations of classical chaos. Dynamical localization 
[3| and phase space localization (scarring phenomena) 
[!,[| have been investigated in those systems. Recently, 
Wilkinson |^| studied the dispersion of the energy in a 
generic non integrable system when the parameter is time 
dependent. He observed that it spreads diffusively and 
found asymptotic theoretical expresions (using random 
matrix theory) for the diffusion constant in the limits of 
large and small velocities of the parameter. Numerical 
experiments performed in a random matrix model were 
in good agreement with those predictions. The physical 
motivation of these studies is the possibility of modelling 
the quantum dissipation of a finite-size system of non- 
interacting fermions, based on that microscopic diffusive 
behaviour of the energy dispersion 

In this paper, we treat the dispersion of the energy in 
a more realistic system: a two dimensional chaotic bil- 



liard with moving boundary. It is not a simple numerical 
task to solve the time dependent Schoedinger equation for 
a 2-D billiard with moving boundary j7j, specially, it is 
very lengthy from the computationally point of view. We 
found a simple shape parameter Hamiltonian expansion 
which enables us to obtain the spectrum (up to second or- 
der) and wave functions of the deformed billiard for defor- 
mations so large as the characteristic wave length in the 
working energy region. Using that expansion and defin- 
ing a shape oscillatory motion, the quantum dynamic of 
a particle inside the billiard reduces to a system of cou- 
pled linear differential equations which can be integrated 
directly. On the other hand, this expansion can be used 
to obtain efficiently static properties of the spectrum like 
curvature distributions or avoided crossing distributions. 

We studied the quantum dynamics of the Bunimovich 
stadium billiard with oscillating walls and calculated the 
diffusion constant D as a function of the boundary veloc- 
ity. Differences with theoretical predictions were found. 
For small velocities of the boundary, the system diffuses 
more than predictions. In this regime, the existence of 
bouncing ball states enables coherent transport of prob- 
ability, so the diffusion is enhanced. On the other hand, 
for large velocities, the system diffuses less than predic- 
tions. We propose phase space localization as the possi- 
ble mechanism to reduce diffusion in this regime. These 
assertions are supported with numerical simulations on 
the stadium billiard where the most localized states in 
the region under study were extracted. In this case, the 
system behaves in good agreement with the theory. 

The outline of the paper is as follows. To make it self- 
contained, in Sec. II we give a short introduction to the 
main results obtained previously for the diffusion of the 
energy in generic chaotic systems. In Sec. Ill we show 
how to expand the Hamiltonian of a 2-D planar billiard 
in powers of the shape parameter. In Sec. IV we give a 
detailed numerical study of the diffusion constant D as a 
function of the velocity of the boundary for the stadium 
billiard. Section IV is devoted to final remarks. We in- 
clude an appendix with the expresions required to obtain 
the parameter Hamiltonian expansion in the stadium bil- 
liard. 
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II. QUANTUM ENERGY DIFFUSION IN 
CHAOTIC SYSTEMS 



Let H(£) be a parametric Hamiltonian related to a 
generic chaotic system for all I. The parameter £ is a time 
dependent function £(t). We restrict to systems with time 
reversal symmetry where the statistical properties of the 
spectrum are well described by the eigenvalues of ensam- 
bles of gaussian orthogonal random matrices (GOE) || . 
At any £{t), the system admits an energy spectrum given 
by the eigenvalue problem for the instantaneous Hamil- 
tonian 



2m 



Mm,r). (2.1) 



The system is prepared at t — in a highly excited eigen- 
state ip v (e(0),r) of the Hamiltonian H(£(0)). The state 
for time t can be expressed in the b asis of eigenstates 
of the instantaneus Hamiltonian (eq. 2A ) , the so called 
adiabatic basis, 

*(r,i) Mm,r) . 



The dispersion of the energy for the considered state is: 
A 2 E(t) = Mt)\\Mm) ~ E(t)) 2 , (2-2) 



where E(t) = (tf = £^ KWIXW*)) is the ex " 
pectation value of the energy as a function of t. 

Wilkinson [6j observed that the dispersion of the en- 
ergy exhibit a diffusive growth when it is averaged over 
many states, that is 

{A 2 E){t) = 2 D t. 

Using random matrix theory for the statistical properties 
of the spectrum, he predicted that the diffusion constant 
D is an universal function of p (the mean energy density) , 
h, £ and the mean value of the off diagonal elements of 
dH/dl that he called a. The dimensionless parameter 
k = p 2 fi£a is a measure of the adiabaticity of the varia- 
tion in the Hamiltonian. He found the following asymp- 
totic form of the diffusion constant for small and large 



D 



5 / 4 7rr(3/4)^-K 3 / 2 if«<l 
if k » 1 



(2.3) 



III. 



HAMILTONIAN EXPANSION FOR 
DEFORMED BILLIARDS 



This section is devoted to obtain a Hamiltonian expan- 
sion for a particle of mass m inside a shape parameter 



dependent planar billiard. For variations 8 = £—£q of the 
parameter around £q, we propose the following expansion 
in the basis of eigcnfunctions at £ , 



H^(£o + 6) ~ (H^ + S 2 H'^/2) 8^ + 5 H'^, 



(3.1) 



with <5 Wi = 1, and 8^ v — Q for ji ^ v. This allows us to ex- 
press the eigenfunctions of the deformed billiard in terms 
of the eigenfunctions at £q. Then, when the parameter 
changes as a function of time, the dynamical evolution 
of the particle can be integrated easily using a standard 
Runge-Kutta method (see Sec. IV). 

Let ( be a smooth close curve defining a planar bil- 
liard. We use a curvilinear coordinate system around the 
boundary with s along £ and z perpendicular to it at s 
(z = on £). Consider now that the boundary is deformed 
and the changes are parametrized by 



r(s,S) = r (s) + z(s,S) n , 



(3.2) 



with i"o(s) the parametric equation for £, and n the out- 
ward normal unit vector to £ at ro(s) (see Fig. 0). In 
the appendix we give z(s, 6) for the stadium billiard. 

Let {^> M (r) ; p — 1,2, ...N} be the eigenfunctions of 
the billiard defined by £ with eigenwave numbers fc M 
around ko (\k^, — ko\ < Ak ~ Perimeter/ Area) . Tak- 
ing h 2 /2m=l, we immediately see from Eq. (|3.l[) that 



H^i^i = k^. The eigenfunctions vanish on £ and can be ex- 
tended outside the billiard in a smooth way; up to second 
order as follows: 



M (r o ( S ) + z n) = z ^t(r (a)) + 0(z 3 ) 
on 



(3.3) 



Moreover, to each function M (r) we associate the scaling 
function cjy^kr/k^) . This family of functions depending 
on the scaling parameter k verifies Helmholtz equation 
with wave number k. 

For small variations of the boundary {5 -C 1), it is valid 
to use perturbation theory to obtain the eigenfunctions 
and eigenvalues of the deformed billiard. With those 
aproximated solutions we will obtain a Hamiltonian ex- 
pansion in powers of 8 which allows us to extend the 
range of deformations to the order of the wave lenght 
(27r/fco). 

Let Vv( r ' ^) b e the eigenfunction of the deformed bil- 
liard obtained from 4>^ (r) by a continuos variation of the 
parameter. As we go to obtain the spectrum up to sec- 
ond order in 8, we express ipf,, in terms of <p„ up to first 
order 



Vv(r,(5) = <^ 



k]j 



with k^{8) the wave number for which ijj^(v,S) van- 
ishes on the deformed boundary (Eq. ( |3.2| )). Then, by 
expanding around £q we have 
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0=^(ro( S ),0) 



r n /k,+z'(s)]^t(v (s)) + 

c P-y <MV"oOO/M , (3.4) 



where primes indicate derivation with respect to £ at £q 
and r n = r.n. In what follows, we omit the argument 
(r (s)) of the normal derivatives and drop terms includ- 
ing derivatives of second order or more. In particular, 
we take </>„(fc M r /fc„) — {k^/k v — 1) d<\> v jdn . Then, 
multipying Eq. ( |3 ,4| ) by <9</> M /3n and integrating on £, 
we have 



2 \l 



(K) 



z'(s) { d -^fds 



(3.5) 



To obtain the last expresion it is necessary to use the 
following quasiorthogonality relation 0] 



dn dn 2kak u ^ 



(kfi k v ) 
(kfj, + k v ) 



0(1) 



Now, multiplying Eq. (3.4) by d(j) v /dn and integrating 
on £, we obtain 



2k v 



z'(s) 



C 



<9n <9n 



ds. 



By perturbation theory the last expresion would be H'] 
however, it is not symmetrical with respect to /i and v. 
Working to first order in 5 it is sufficient to replace the 
factor [ky, + fc„)/2fc„ by one (we assume that \k^ — k u \ <C 
kg); but to second order it is not the case. Consider 
two states at the same distance from k^; that is k± — 
k^ ± e. Then, the factor in question (to first order in 
e) is 1 =F e/2fc M . To solve the problem, we define the 
symmetrical factors 



A. 



(fc P + K)/2k 



which give the same asymmetrical contribution to the 
state n from states ±. On the other hand, we stress that 
for generic non-integrable billiards and for any non trivial 
deformation (a dilation is trivial) the interaction is of long 
range. So, to have into acount the finite dimension of the 



basis we suggest to multiply the elements of Eq. (3.5) by 
the following cut-off: 



exp -2 (ki - kty/(k Ak) 



Finally, non diagonal Hamiltonian elements acquire the 
expresion 



-Cf^ Apu (p z'{s) 



dn dn 



ds. 



Working as before, to obtain H'^ we multilply equa- 
tion = d 2 ijj fl /d£ 2 (r (s), 0) by d<j)^/dn and integrate 
it on £. The result is 



z"(s)(^fds 



Figure g compares the approximated spectrum obtained 
from eq. (3.1) with the exact one for the Bunimovich 
stadium billiard; the agreement is excellent. 



IV. NUMERICAL RESULTS 

We have consider the desymmetrized stadium billiard 
with radius r and straight line of lenght a. The bound- 
ary only depends on the shape parameter £ — a/r (the 
area is fixed to the value 1 + ir/4). The parameter £ os- 
cillates harmonically around £ = 1 with frequency uj and 
amplitude a, that is, 

£{t) = 1+ a sin(wt). 

Sec. II shows universal expresions for D in the asymp- 
totic limits of k; therefore, n is a good quantity to char- 
acterize the quantum dynamics of a specific chaotic sys- 
tem. In order to calculate the coefficient D, we solved 
the time dependent Schroedinger equation for different 
initial conditions 



ih 



dt 



H(£(t)) 



(4.1) 



A very efficient way (computationally) to solve 4.1 is us- 
ing the parameter expansion of the Hamiltonian obtained 
in Sec. III. The solution ^(r, t) is expanded in the basis 
of eigenf unctions at £(Q) = 1, 



*(r,t) = <Mr). 



(4.2) 



Although ^(r) satifies boundary conditions for £(0) = 1 , 
it does not cancel outside the boundary (see Eq. |3.3| ). 
Of course, the amplitude of the oscillations is limited to 
values of the order of the wave lenght. After replacing 
4.2 in 4.1, we obtain the system of differential equations 



This system was integrated using a standard fourth- 
order Runge-Kutta method, and the dispersion of the 
energy is evaluated using (2.2) at times ( = 2 n ir/w 
(n = 0,1,2,...). We take h = 1 and m = 1/2. Then, 
we used a basis of 103 odd-odd eigenfunctions around 
fco = 47.3 at £ = 1 (see Fig. ||). Time evolutions were 
stopped before the wave functions had spread into the 
tails of the spectrum. 

As an example, Fig. |^ shows the dispersion of the 
energy < A 2 E > for u — 0.4 and a = 0.05. The av- 
erage is over ten initial states. They were chosen near 
the center of the spectrum. It is observed clearly a linear 
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spreading (diffusive behavior) of the dispersion of the en- 
ergy for the first oscillations of the boundary (the slope 
being equal to 2 D); later, the spread goes slowly and 
eventually saturates. 

Figure ^ shows a log-log plot of the difussion constant 
D as a function of n which is proporsional to the mean 
value of the boundary velocity < \£\ >— 2aLu/n. Simula- 
tions were performed for uj = 0.05, 0.1, 0.2, 0.4, 1, 2, 5, 
10, 20 and 40, with a = 0.05 (filled circles). In all the 
cases the average was over ten initial conditions. Er- 
rorbars were calculated using the standard deviation of 
the average in the ensamble of initial states. The figure 
shows theoretical predictions in solid straight lines for a 
gen eric chaotic system with time reversal symmetry (Eq. 
2.3), and we can observe significant differences with the 
numerical data. 

In order to directly see the influence of localized eigen- 
functions on the quantum dynamics of the system, we ex- 
tracted from the spectrum the most localized states. This 
procedure is very simple by using the Hamiltonian expan- 
sion presented in Sec. III. In order to extract state j, we 
remove simply row and column j from the Hamiltonian 
matrix. If at I = 1 this state collide in an avoided crossing 
it is necessary to do first the transformation explained in 
reference Q . We have removed bouncing ball states with 
k = 46.4589, 47.1943, 47.4220, and 47.7548, all of them 
at £ = 1 (see Fig. ||). The resulting spectrum is showed 
m Fig. |. In the same way as before, we have calculated 
the difusion constant for u> — 0.2, 0.4, 1, 5, 10 and 20 
with a = 0.05 , and for uj = 0.2, 0.4, 10, 20, 40 and 60 
with a = 0.09. In this case, the results are very close to 
theoretical predictions except for the highest values of uj 
(see Fig. |). 



V. FINAL REMARKS 

We have studied the spreading of the energy dispersion 
in the Bunimovich stadium billiard with oscillating walls. 
To the best of our knowledge this is the first computation 
of diffusion on a quantum chaotic system. We arrived to 
that possibility by using a Hamiltonian expansion which 
we develop in this work. 

Numerical simulations show that the energy spreads 
diffusively for a number of oscillations (linear behaviour) 
and then saturates. This saturation phenomenon is asso- 
ciated to the fact that the eigenfunctions of the evolution 
operator for one period (Floquet states), which are local- 
ized in energy, inhibit the spread of energy for long times 
jll],[l2| . We calculated the diffusion coefficient (the slope 
of the initial linear behaviour) for different velocities of 
the boundary, observing significant differences with the- 
oretical predictions based on random matrix theory. 

In the adiabatic regime (slow boundary velocities), the 
system diffuse more than predictions. We associate these 
differences to the existence of states strongly localized 



in phase space (their Wigner or Hussimi distributions 
are strongly localized). In particular, we refer to bounc- 
ing ball states which are strongly localized in momentum 
space. The interaction between these states and generic 
chaotic ones is much smaller than the typical interac- 
tion between generic states. This fact is reflected on 
the occurrence in the spectrum of structures surviving 
parametric variation; that is, straight lines interrupted 
by very small avoided crossings (see Fig. |^). Then, 
having in mind that in this regime Landau-Zener transi- 
tions p3P] are the mechanism of diffusion, those straight 
lines behave like channels allowing a coherent transport 
of probability to a wide range of energy in each oscilla- 
tion. In order to confirm numerically the above explana- 
tion we have extracted from the spectrum those bouncing 
ball states living in the region under study. In this case, 
the system diffuse according to the theory. Certainly, 
classical billiards with a bouncing ball continuous family 
of orbits display anomalous diffusion too Hj. 

As the velocity increases, Landau-Zener probability 
transitions across small avoided crossings approach ex- 
ponentially one. Then, if the system starts in a bouncing 
ball state, it moves essentially up and down by a straight 
line without mixing. And this situation mimics the evo- 
lution of a regular system where diffusion is known to be 
very small fill . In the spectrum without bouncing balls 
states, observed diffusion agrees with theory (in Fig. || 
the solid line with highest slope) reasonably well. How- 
ever, for high velocities simulations and theory fall apart 
again showing the beginning of saturation for the diffu- 
sion process. 

Recently we have proposed Q the elimination of 
avoided crossings as the natural mechanism to uncover 
localized structures embedded in the eigenfunctions of 
chaotic Hamiltonian systems. Precisely, we have shown 
that many scars of short periodic orbits are uncover in 
the stadium billiard transforming the parametric spec- 
trum in a set of smooth curves which cross among them. 
As we mention above, the elimination of avoided cross- 
ings can be carried out dynamically by increasing the 
boundary velocity; that is, when the probability transi- 
tion is practically one. In the light of this, saturation of 
diffusion is expected in chaotic systems for high velocities 
as a consequence of these localized structures associated 
to short periodic orbits. 

In conclusion, we have verified that asymptotic ex- 
presions (2.3) for the energy diffusion in a chaotic sys- 
tem work very well in a wide range of boundary veloci- 
ties except for high velocities where diffusion saturates. 
Moreover, if the system includes a fraction of localized 
states like bouncing ball (or states living in a classically 
regular region which are connected to chaotic states by 
tunnelling), diffusion will be anomalous in the adiabatic 
regime. 
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APPENDIX A: 

In this appendix we show the parametrization used to 
describe the deformation of the desymmetrized stadium 
billiard. The area is fixed to the value A = l + ir/4, so the 
boundary only depends on the shape parameter t = a/r. 
A point ro = (x,y) on the boundary at £q = 1 is given 
in terms of the curvilinear coordinate system defined in 
Sec. II (see Fig. @) 

[ > " \ l + sin(s- 1) if s > 1 . 
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y(s) 



1 if s < 1 

cos(s — 1) if s > 1 



Then , if the deformed boundaries are described by Eq. 
( |3.2|) , it is a geometrical problem to show that z(s, S) ~ 
5 z'{s)+5 2 z"{s)/2, with 



z'(s) = 



i 

' 2 A 



if s < 1 



(1 _ i )sin(s-l) - i if s > 1 



and 



3 

O 5 



if s < 1 



z"{s) 



i-l + (l-^) 2 sin 2 (. S -l)+ 
, 2^-i(!- Ta) sin(s-l) if s > 1 
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FIG. 1. Schematic figure showing the curvilinear coordinate system on the boundary of the stadium billiard. In dashed line 
it is shown a deformation of the billiard. The variations of the boundary are described by the function z(s, 8). 
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FIG. 2. Approximated spectrum (solid lines) obtained from eq. (3.1) is compared with the exact one (dots) for the stadium 
billiar with fixed area. 
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FIG. 3. Dispersion of the energy (A 2 E) vs t (time) for ui = 0.4. The average is over ten initial states. The solid line is the 
best linear fit used to estimate the diffusion constant (the slope corresponds to 2D). 
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FIG. 4. Log-log plot of D vs. k. The (•) correspond to the stadium billiar for uj = 0.05, 0.1, 0.2, 0.4, 1, 2, 5, 10, 20 
and 40. For the case in which bouncing ball states (see text) were removed, the (o) correspond to u> — 0.2, 0.4, 1, 5, 10 and 
20 with a — 0.05 and, the (o) correspond to ui = 0.2, 0.4, 10, 20, 40 and 60 with a = 0.09. It is showed in solid lines the 
asymptotic theoretical prediccions for a generic chaotic system with time reversal symmetry (eq. 2.5). Errors smaler than the 
symbols are not ploted. 
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